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Abstract 

Measurements of the drift velocities of electrons and holes as functions of electric field and tem¬ 
perature in high-purity n- and p-type silicon with <100> lattice orientation are presented. The 
measurements cover electric field values between 2.5 and 50 kV/cm and temperatures between 233 
and 333 K. For both electrons and holes differences of more than 15% are found between our 
<100> results and the <111> drift velocities from literature, which are frequently also used for 
simulating <100> sensors. For electrons, the <100> results agree with previous <100> measure¬ 
ments; however, for holes differences between 5 and 15% are observed for fields above lOkV/cm. 
Combining our results with published data of low-field mobilities, we derive parametrizations of 
the drift velocities in high-ohmic <100> silicon for electrons and holes for fields up to 50kV/cm, 
and temperatures between 233 and 333 K. In addition, new parametrizations for the drift veloci¬ 
ties of electrons and holes are introduced, which provide somewhat better descriptions of existing 
data for <111> silicon than the standard parametrization. 
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1. Introduction 

Simulations of silicon sensors and the extraction of parameters from measurements, rely on 
the knowledge of a number of material parameters, for example the drift velocities of electrons 
and holes as functions of electric field and temperature in high-resistivity silicon. Presently, the 
information on drift velocities for silicon with <100> lattice orientation is quite limited. However, 
<100> silicon is widely used in high energy physics detectors as the defect concentration at the 
surface is small compared to <111> silicon. Therefore, results from <111> silicon are typically 
also used for the simulation and analysis of <100> data. We have investigated how well we can 
simulate the measured current transients in non-irradiated pad sensors built on <100> silicon 
and found that by using the data of Refs. mm we were unable to describe the measurements. 
Therefore, we determined the drift velocities of electrons and holes in silicon of <100> lattice 
orientation as functions of the electric field and the temperature using two different methods. 


2. Sensors investigated and experimental setup 
2.1. Sensors investigated 

Three different sensors were investigated, p^nn+ and n+pp+ pad diodes, two manufactured by 
Hamamatsu (HPK) j3] and one by CiS [4], on <100> float-zone silicon. Physical parameters of 
the diodes are given in Table The doping concentration was obtained from capacitance-voltage 
measurements. The mean oxygen concentration of the diodes was 5 — 10 • 10^® cm“^. The active 
thicknesses of the diodes were obtained from the mechanical thickness, which was measured using 
a caliper and dielectric measurements, minus the depths of the implants. For the HPK diodes 
the depths of the n+ and p+ implants, determined using spreading resistance measurements by 
SGS Institut Fresenius GmbH |5], were used. As we have no measurements of the implant depths 
for the GiS diode, the values of the HPK diodes were used. The values determined from the 
capacitance above full depletion were used as a cross check. The error of the active thickness is 
estimated to be 1 %. To allow injection of light from both sides of the pad diodes, the junction 
(front) side contact is a square aluminum ring, which leaves free the center of the diode, and the 
opposite (rear) side contact is an aluminum grid. 


Vendor 

Type 

w 

[pm] 

Udep 

[V] 

^eff 

[cm“^] 

A 

[mm^j 

Cfd 

[PF] 

{^Tnin} 

[kV/cm] 

i^max') 

[kV / cm] 

HPK 

n 

200 ±2 

87.5 ±2.6 

2.9- 10^^ 

4.4 

2.7 

5.0 

50.0 

HPK 

p 

200 ±2 

115.0 ±3.5 

3.8-10^^ 

4.4 

2.7 

6.0 

50.0 

GiS 

n 

287 ±3 

48.2 ± 1.4 

0.8- 10^^ 

24.5 

9.2 

2.4 

34.8 


Table 1: Physical properties of the investigated pad diodes; w denotes the active thickness, Udep 
the depletion voltage, Ngff the effective bulk doping, A the pad area, Cfd the capacitance above 
full depletion, and {Emin) and (Emax) are the minimum and maximum mean electric fields in the 
diodes where the drift velocity was measured. 


2.2. Experimental setup 

In order to determine the drift velocities of electrons and holes separately the transient current 
technique (TGT) was employed, with charge carriers produced by pico-second laser light with 
different wavelengths. Most of the measurements were performed at the Hamburg setup, which 
is described in more detail in Refs. 017]. Additional control measurements were performed at 
the TGT setup at the GERN-SSD Lab |H|. Gurrent transients generated by laser pulses with a 
EWHM = 50 ps and a pulse frequency of 200 Hz were recorded via 3 m RG58 cable, a bias-T, 
and an amplifier 0 using a Tektronix DPO 4104 oscilloscope with a bandwidth of 1 GHz and 
a sampling rate of 5GS/s. To reduce the effects of noise 512 transients were averaged at each 
measurement step. 
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Two lasers [TO] with different wavelengths were used: a red laser with light of A = 675 nm 
wavelength for front- and rear-side illumination, and an infrared laser with A = 1063 nm for 
front-side illumination. The absorption length of 675 nm light in silicon is A « 3.3 pm at room 
temperature (RT), and charge carriers are created only close to surface of the diode. Thus only 
one carrier type drifts through the whole diode and generates most of the transient current. For 
illumination of the p+n junction of a p+nn^ diode, the generated holes are collected immediately 
at the close-by p^ contact, while the electrons drift through the n bulk to the n+ contact. Likewise, 
for the illumination of the n+ contact the electrons are collected immediately, while the holes drift 
through the bulk to the p+ contact. The absorption length of the infrared laser light of A « 1 mm 
at RT is much larger than the active thickness of the investigated diodes. Therefore, the electron- 
hole pairs are generated approximately uniformly along the path of the light and both electron 
and holes contribute equally to the current transient. 

The measurements were performed at different temperatures for bias voltages of Udep < t4ias ^ 
1000 V in steps of 10 V. For the HPK n-type diode, the measurements were performed at seven 
temperatures: T = (233, 253, 263, 273, 293, 313, 333) K. For every value of Utias and T, the 
measurements were performed for all three laser configurations: front illumination with the red 
laser, rear illumination with the red laser, and front illumination with the infrared laser. 


3. Analysis methods 

3.1. Transient current simulation 

A simulation of charge transport including the diffusion of charge carriers was used to describe 
the TCT measurements. A more detailed description can be found in Refs. [HIT]. The following 
description of the simulation is partly taken from Ref. m. Electron-hole pairs are generated on 
a grid with Ax = 50 nm spacing according to an exponential function with the light-attenuation 
length A(T). The temperature dependence of A(T) is modeled according to Ref. [T^]. The charge 
carriers are drifted in time steps of At = 10 ps, taking into account diffusion by a Gaussian with 
variance = ^ ■ fif,- ks ■ T/qq) ■ At for electrons, and similarly for holes. The Boltzmann 
constant is ks, the absolute temperature T, the elementary charge qq, and the field-dependent 
electron and hole mobilities fj-e and The current induced in the electrodes in the time interval 
between i ■ At and (i -|- 1) • At is calculated according to J®*'" = qo/At- J2j [i^i+i j ~ ^tj) ~ 
j — ) ], where Nfj is the number of electrons and N^j the number of holes at the grid 

point j at time i ■ At. Effects like charge trapping or charge multiplication are not considered. 
A constant space charge density go ■ N^ff is assumed, resulting in a linear electric field. It has 
been verified, that the current transients are hardly affected if the electric field simulated with 
SYNOPSYS-TCAD [T3] for realistic doping distributions at the implants is used instead. 

The transfer function R of the electronics circuit is determined using the convolution theorem 
J-{f ®g} = iF{f} ■ IF{g} of Fourier transforms. The measured TCT pulse has a sampling interval 
of 200 ps, and a spline interpolation of the measurement is used to obtain values for the 
10ps time steps of the simulation. The transfer function is obtained by i? = ^ ■ 

The convolved simulated signal is then obtained by S'!*™ = compared to the 

measurements. The method is described in more detail in Ref. m- 

3.2. Mobility parametrizations 

The Caughey-Thomas (CT) parametrization |14[ [T] is generally used to describe the field 
dependent mobility in high-ohmic silicon 



with the low field mobility p-o, the saturation velocity Vsat, and a phenomenological parameter /3. 
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When comparing different mobility measurements reported in the literature, we noticed that 
for high electric fields l/fj,{E) is to a good approximation a linear function of electric field E, with 
a very sudden transition from the constant low-field mobility fj,o to the linear E dependence. For 
illustration Fig.shows l/fi{E) for electrons and holes for <111> silicon at T = 245 and 300 K 
obtained by digitizing figures of Ref. m, and a fit by the new parametrization (KS) presented in 
this paper 




E<Eo 

+ V-S -{E-Eo) E>Eo. 


( 2 ) 


At E = Eq the inverse mobility changes from a constant value to a linear increase. In order to 
have a function with a continuous derivative, one may use the function l/fi{E) = l/jiQ + (E — 
Eq) ■ (1 -I- tanh(£' — EQ))/{2vsat), with the units of V/m for E. The maximum deviation of this 
equation from Eq. |^is 10“^%. We note that for Eq = 0 the new parametrization has the same 
field dependence as the Trofimenkoff (Tr) parametrization |16| : 

Tr 

E^^{E) = ( 3 ) 

The physics motivation for the parametrizations is: at low electric fields the drift velocity Vd is 
small compared to the thermal velocity Vth- Assuming an effective scattering rate 1/t lattice ■, the 
drift velocity Vd{E) = {{qq ■ E/m*) ■ dt = (qq ■ nattice/m*) ■ E, and the low-field 

mobility is constant: hq = qq ■ Tiatuce/ tti*, with the effective mass of the charge carriers to*. At 
high electric fields the charge carriers can acquire an energy Cemission, which is sufficient for phonon 
emission. The phonon emission rate 1/TemissioniE), which increases with increasing electric field, 
is responsible for the additional term cx if in Eq. leading to saturation of the drift velocity. We 
introduce a threshold field Eq for the phonon emission. 
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Figure 1: Inverse of the mobilities for electrons and holes as functions of the electric field for 
<111> silicon at T = 245 and 300 K. The data have been obtained by digitizing the points of 
Figs. 5 and 7 of Ref. [T5]. The curves are fits to the data using Eq. 


We refer to the Appendix 1 for a comparison of the fit results for the two mobility parametriza¬ 
tions to the electron and hole data for <111> silicon of Ref. m- Here we only mention that at 
temperatures above 245 K the new parametrization provides a somewhat better description of the 
data. At lower temperatures none of the parametrizations provides a good description. For the 
quantitative comparison we use the root-mean-squared relative deviation 
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As an example, for the data at T = 300 K the values of rdev using the CT parametrization 
are 0.41% for electrons and 0.42% for holes, whereas they are 0.17% and 0.28% for the KS 
parametrization. 

For the hole mobility the quality of the description of our data could be improved significantly 
by adding a quadratic term to the parametrization of the inverse mobility of Eq. 


^/Th{E) 


VmS E <Eo 

i/^^^+b.{E-Eo)+c-iE-Eof E>Eo. 


( 5 ) 


The parameters of the KS mobility model were determined for the temperatures at which the 
data were taken, as well as for the entire temperature range. For the parameter c no temperature 
dependence has been assumed. The temperature dependencies of the parameters was modeled 
similar to Refs. [manz] as a power law 

/ T [Kl \ 

pari{T) = par,(T = 300K) • , (6) 

with pari(T = 300K) denoting the parameter values at 300K and the power 


3.3. Fit method 

For determining the drift velocities of electrons and holes the quantity 

n 

x" = const (7) 

i=l 

was minimized, where are the measured current values, and the simulated ones after 

convolution with the transfer function R. As discussed in Section]^ fits have been made for 
different data sets. For each sensor the transfer function R, determined from the measurements 
taken with the infrared laser at T = 313 K and Uuas = 1000 V, was used. For determining R the 
drift velocities for electrons and holes have to be known. Therefore, R had to be determined in an 
iterative way: for a given set of measurements R has been determined using the Fourier transform 
method described in Ref. m- Using this transfer function a set of mobility parameters has been 
obtained by a fit of the entire data set for a given sensor, meaning all transients at bias voltages 
above depletion at all temperatures for all three laser configurations. This procedure was repeated 
until stable results have been reached after about 3 iterations. It has been checked that the final 
results do not depend on the starting values of the mobility parameters. 

As shown in more detail in Ref. m, the entire data set is well described by the fits. To 
illustrate the quality of the description of the data, we show four examples in Fig.[^ the transients 
for Ubias = 1000 and 150 V at T = 313 K for the front- and the rear-side illumination with the 
675 nm laser light. Details of the pulse shape, including signal reflections at connectors and the 
amplifier, are well described. 


3 . 4 . Time-of-flight method 

The well established time-of-fiight (tof) method, used e.g. in Ref. [15], was used as a cross¬ 
check of the results of the fit by the simulated transients discussed above. As shown below, the 
tof method can be extended to the situation of non-uniform electric fields, if the inverse mobility 
is a linear function of the electric field, as it is the case for the parametrization of Eq.|^ below 
and above the field Eq. The data with front-side injection of the 675 nm laser were used for 
determining the electron drift velocity, and the ones with rear-side 675 nm light injection for the 
hole drift velocity. For the time-of-fiight, tto/j the time difference between the maximum slopes at 
the rise and at the fall of the pulse has been used. Using the simulation discussed above, we have 
verified, that this definition of Lof agrees with the calculated transit time {p{E) • E{x)) Ax 
to within 100 ps for Ubias > 150 V. 
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(c) (d) 


Figure 2: Comparison of the measured current transients (crosses) with the simulated ones (solid 
lines) for the HPK n-type diode using 675 nm laser light measured at T = 313 K. (a) Front-side 
illumination at Ubias = 1000 V (electron signal), (b) rear-side illumination at Ubias = 1000 V (hole 
signal), (c) front-side illumination at Ubias = 150 V (electron signal), and (d) rear-side illumination 
at Ubias = 150V (hole signal). 


As the electric field, E{x), depends on the position in the sensor, x, also the drift velocity 
Vd{x) depends on x. The electric field, Etof, corresponding to the drift velocity Vd{Etof) = w/ttof 
obtained by the time-of-flight measurement, has to be determined. As shown below, E^af = 
{\/E[x))~^, with () denoting the average over a;, if the inverse of the mobility is constant or linear 
in E. These conditions are satisfied if the mobility can be described by Eq. [^and the electric field 
in the entire sensor is either smaller or larger than Eq. 

Here the argument: from l/fi{E) = ai + a 2 ■ E{x) follows l/vd{E) = ai/E{x) + 02 , and from 
= jQ{l/vd{E))dx = {ai/E{x) + a 2 )dx = w ■ {ai{l/E{x)) + a 2 ) and Uof = w/vd{Etof), 
follows 1/Etof = l/{E{x)). For uniform doping Etof = {Emax - Emin) / '^'<^{Emax/E^in) , where 
Emax and Emin denote the maximum and minimal fields in the sensor, respectively. 

For a sensor of 200 pm thickness and 90 V depletion voltage, for bias voltages of 100, 150, 250, 
500, and 1000 V the values for {E) are: 5, 7.5, 12.5, 25, and 50kV/cm. The corresponding values 
for {1/E)~^ are 3.06, 6.49, 11.9, 24.7, and 49.9kV/cm. The differences between (E) and {1/E)~^, 
in particular for bias voltages close to the depletion voltage, are significant. 

As l/vd{E) = ai/E{x) + a 2 is only an approximation for Eq.[^ we also compared the measured 
values of ttof with the values of ttof = Jq (/^(T?) • E{x)) ^dcc using the values for fJ,{E) from the 
the parametrization Eq.j^for electrons and Eq. [^for holes. However, the differences to the method 
described above are small, and the linear dependence of the inverse mobility provides an adequate 
approximation for the tof method even close to the depletion voltage, where the variations of E{x) 
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and Vd{x) are large. 


4. Results 

4-.1. Fit and time-of-flight results 

We first fitted separately the data of the HPK n-type diodes for the individual temperatures 
using the parametrization of Eq.j^for electrons and of Eq.j^for holes. As the minimum investigated 
mean electric field was 5 kV/cm, our data are not sensitive to the constant low-field mobility of 
electrons below 1.8kV/cm at room temperature. For reaching lower field values thicker sensors 
with higher resistivity as used by Canali et al. are required. The results are given in Table[^ 



Electrons 

Holes 

r[K] 

/Jo [cm^/Vs] 

Vsat [10 ''cm/s] 

/J .0 [cm^/Vs] 

& [10 '^s/cmj 

c[10-^®s/V] 

Eq [V/cm] 

233 

2333 

1.127 

842.6 

1.040 

-4.899 

2096 

253 

2013 

1.103 

697.9 

1.024 

-4.608 

2123 

263 

1882 

1.089 

631.9 

1.019 

-4.558 

2379 

273 

1753 

1.078 

576.8 

1.003 

-4.293 

2476 

293 

1523 

1.055 

486.1 

0.9705 

-3.728 

2667 

313 

1327 

1.031 

408.8 

0.9390 

-3.229 

3313 

333 

1157 

1.014 

349.2 

0.8938 

-2.443 

3980 


Table 2: Results of the fit of the simulated to the measured transients at the individual tempera¬ 
tures T for the HPK n-type sensor, using Eq.j^for the parametrization of the field dependence of 
the mobility of electrons, and Eq.l^for the one of holes. 


We then performed a global fit to the entire data set: the HPK n-type sensor illuminated with 
red laser light from the front and rear side as well as with the infrared laser, for temperatures 
between 233 and 333 K and voltages between 100 and 1000 V in steps of 10 V. The corresponding 
mean electric fields are between 5 and 50 kV/cm. Overall, more than 10® data points, each 
one consisting of 512 measurements, were used for the fit. For every iteration the transient for 
every data set had to be simulated. For the temperature dependence of the parameters, the 
parametrization of Eq.j^has been used. The results are given in Table 

The entire data set is well described by the global fit. To illustrate the quality of the de¬ 
scription, Fig.shows the transients for different bias voltages for the lowest (233 K) and highest 
(333 K) temperatures for front- and rear-side-illumination with the 675 nm laser light. The biggest 
differences between fit and data are observed in Fig.|3b| around 10 ns at a bias voltage of 100 V. 
As this voltage is only about 10 V above the depletion voltage, the charge carriers are generated 
in a low field region, where the plasma effect |18l [19] may cause an additional delay and diffusion 
for the holes. 

In Fig. 1^ we compare for 233 and 333 K the field dependencies of the electron and hole ve¬ 
locities from the fit at this particular temperature to the ones from the global fit and from the 
tof measurements. The drift velocities are shown on top, the ratio of the velocities from the tof 
measurement to the fit at the particular temperature in the middle, and the ratio of the velocities 
from the global fit to the fit at the particular temperature at the bottom. We note that, with 
the exception of electric fields below 5 kV/cm where the tof measurements become inaccurate, the 
results from the tof measurement and the fit at a single temperature agree within about 1 %. This 
gives an idea of the accuracy of the analysis methods. The comparison of the global-fit results to 
the single-temperature-fit results shows that, for the electric fields where the measurements were 
made, the agreement is within about 2 %. However, extrapolating to lower and higher field values, 
the differences, in particular for holes at T = 233 K, rapidly increase. 
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pari{T = 300 K) 


Electrons 


1430 cm^/Vs 

-1.99 

Vtat 

1.05 • 10^ cm/s 

-0.302 

Holes 


457 cm^/Vs 

-2.80 

b 

9.57- 10-« s/cm 

-0.155 

c 

-3.24- 10"^^ s/V 

- 

Eo 

2970 V/cm 

5.63 


Table 3: Parameters for the mobility for <100> silicon obtained from the fit to the data of the 
HPK n-type diode using Eq. [^for the electron mobility, Eq. [^for the hole mobility, and Eq. |^for 
the temperature dependence. The mean electric field values range from 5 to 50kV/cm, and the 
temperatures from 233 to 333 K. The uncertainty is estimated to be 2.5% for the electron and 
5 % for the hole mobility. 






(c) (d) 

Figure 3: Comparison of the measured current transients (crosses) with the simulated ones (solid 
lines) for the global fit discussed in the text for the HPK n-type diode using 675 nm laser light and 
bias voltages between 100 and 1000 V. (a) Front-side illumination at 233K (electron signal), (b) 
rear-side illumination at 233K (hole signal), (c) front-side illumination at 333K (electron signal), 
and (d) rear-side illumination at 333K (hole signal). 


4-.2. Comparison of different diodes 

In order to check the results from the HPK n-type diode, we also performed measurements 
on a 290 pm-thick n-type diode produced by CiS and a 200 pm-thick p-type diode produced by 
HPK, with the parameters shown in Table[2 We expect different systematic errors for the different 
diodes. As an example, the full-width-at-half-maximum of the transfer function is about 0.6 ns for 
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drift velocity @ T = 233 K drift velocity @ T = 333 K 



Figure 4: Top: Comparison of the drift velocities of electrons (black) and holes (red) from the tof 
measurements (crosses), the fits at a single temperature (solid lines), and the global fit (dashed 
lines). Middle: Ratio of the tof to the single-temperature-fit results, and bottom: Ratio of the 
global-fit to the single-temperature-fit results, (a) Results for 233 K, and (b) for 333 K. 


the HPK diodes, and about 1.0 ns for the CiS diode, which has a 5.6 times larger area and is 40 % 
thicker. 

As an illustration we show the comparisons of the drift velocities for electrons and holes as 
function of the electric field at 253, 273 and 293/313 K. Fig. 5a shows the ratios of the drift 
velocities for the CIS n-type diode relative to the HPK n-type diode. For both electrons and holes 
the results agree within 1% for fields above 2.5kV/cm, where our measurements are sensitive. 
Fig. 1^ shows the comparison of the HPK p-type to the HPK n-type diode. The electron velocity 
is systematically higher by about 2 — 3 % and the hole velocity higher by about 3 % at a voltage 
of 6kV/cm, and lower by about 1 % at 50kV/cm. The measurements are not sensitive to fields 
below 2.5kV/cm, where the differences are larger. We conclude that over the entire electric field 
range of our measurements, the drift velocities agree to better than 4 %. 

We note that we did not achieve a good description of the transients of the HPK p-type diode 
with the electric field used so far: a linear function from the diode surface at the pn-junction 
side to the other surface. However, a good description was achieved by introducing field-free 
regions of 1 to 2 |a.m at the diode surfaces. In the simulation charge carriers, which are produced 
in the field-free regions, reach the electrodes or the field region by diffusion and thus change the 
transients. For more details we refer to Appendix 2. We also note that the introduction of the 
field-free regions hardly affects the values obtained for the drift velocities from the fits. 

From the differences between the two methods and between the fit results, we estimate for 
fields between 2.5 and 50kV/cm an uncertainty of 2.5% for the electron and of 5 % for hole drift 
velocities for the values of the parameters given in table 


4-3. Comparison with published data 

We compare our results for <100> silicon, presented in Table to the results of Jacoboni et 
al. [T] for <111>, of Becker et al. [5] for <100>, and the TOF results of Canali et al. [TS] for 
<100> silicon. We also compare to the Jacoboni <111> parametrization [T], because it is widely 
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Figure 5: Ratios of the measured drift velocities of electrons and holes at temperatures of 253, 
273 and 293/313 K for (a) CIS n-type diode to HPK n-type diode, and (b) HPK p-type diode to 
HPK n-type diode. The vertical lines denote the maximum and minimum mean electric fields in 
the diodes. 


drift velocity @ T = 245 K drift veiocity @ T = 300 K 



Figure 6: Comparison of the drift velocities from the global fit results of data from the HPK 
n-type diode to published data. On the top the drift velocities and on the bottom their ratios are 
shown for (a) 245 K, and (b) 300 K. 


used for <100> silicon due to the lack of publications on <100> drift velocities. Figure shows 
our results and the published values at 245 and 300 K. On top the drift velocities are shown, and 
at the bottom the ratios to our global-fit results from the HPK n-type diode. 

At both temperatures the Jacoboni <111> drift velocities for electrons are 10 — 20% higher 
than our <100> results. At 245 K the Jacoboni <111> drift velocity for holes is 2 — 10% lower 
than our <100> results. At 300 K the differences, in particular for electric fields below lOkV/cm, 
are smaller. 

When comparing our results to the values published by Becker et ah, one notices that the 
electron drift velocities are very similar, whilst the drift velocities for holes are ~ 15 % below 
our results at high electric fields. However, the measurements of Becker et al. were made at 
(E) = (3.8, 7.1, 17.9) kV/cm, where the differences are less than 6%, and the large differences are 
outside of the region of their measurements. 

When comparing the Canali et al. tof values to our fit results, one notices that the electron 
drift velocities are up to 20 — 25 % higher for fields below 2 kV/cm. However, these field values 
are outside of the range of our measurements. The hole drift velocities determined by Canali et 
al. are similar to our results at low electric fields, but at high electric fields they are up to 15 % 
lower. 
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pari{T = 300K) 


Electrons 

Mo.Jac 

1530 cm^/Vs 

-2.42 

<at 

1.03 • 10' cm/s 

-0.226 

Holes 

f^0,Jac 

464 cm^/Vs 

-2.20 

b 

9.57- 10-® s/cm 

-0.101 

c 

-3.31 • lO-^'^s/V 

- 

Eo 

2640 V/cm 

0.526 


Table 4: The mobility parameters for <100> silicon of Eq. for the electron mobility, Eq. 
for the hole mobility, and Eq. j^for the temperature dependence. The values of Atg’jac 
corresponding values have been taken from |T], and were fixed in the fit to the measurements 
of this manuscript in order to describe the mobility for electric field values up to 50 kV/cm and 
temperatures between 233 and 333 K. However, for E > 2.5 kV/cm we recommend using the 
parameters of Table [3| 


/./. Parametrization of the drift velocities for fields up to 50 kV/cm 

In order to provide a parametrization of the drift velocities of holes and electrons in <100> 
silicon, which can be directly used in simulations, we have combined the low-field data for <111> 
silicon of Ref. |T] with our data. The low-field mobility in silicon is isotropic as shown in Ref. [T]. 
For electrons the parametrization of Eq. for holes the one of Eq. and for the temperature 
dependencies the one of Eq. [^have been used. For the fit the values of po and at of Ref. [1] were 
fixed. The results are given in Table 

We want to stress that we did not observe the expected constant part of the electron mobility. 
The transition region of the electron mobility from a constant value to a linear increase, corre¬ 
sponding to velocity saturation, can not be described by the simple model presented in Eq. 
But, at least one additional parameter similar to the quadratic term c for the hole mobility, Eq. 
is needed for an adequate description. Unfortunately, our data in the transition region between 
about 0.5 and 2.5kV/cm is not sufficient to determine such a parameter precisely. 

For fields between 2.5 and 5kV/cm differences of up to 9% to the fit without fixing the low- 
field mobility are found. There is no error estimation given for the values of Ref. [T]. However, 
the authors of Ref. [1] note the error to be comparable to 5 %. Therefore, we recommend using 
the parameters of Table for E > 2.5kV/cm. 


5. Summary 

The current transients in p- and n-type-silicon pad diodes from electron-hole pairs produced by 
sub-nanosecond laser light with a wavelength of 675 nm injected from both sides, and by laser light 
of 1063 nm injected from the junction side, were measured. From these data the drift velocities of 
electrons and holes in high-ohmic silicon with <100> lattice orientation for electric fields between 
2.5 and 50kV/cm and temperatures between 233 and 333 K have been determined using two 
methods: a fit of simulated transients to the measured transients, and a time-of-flight method. 
The results of both methods agree within 2 %. New parametrizations of the field dependencies of 
the drift velocities have been used, which also provide a somewhat better description than previous 
parametrizations of published results on drift velocities in <111> silicon. The measurements were 
performed in a range of intermediate temperatures which are of actual interest for silicon detectors 
in high energy physics. 

For electrons our results are similar to published results. For holes at high electric fields, 
significant differences to published results are found. The differences to the parametrization for 
<111> silicon of Jacoboni et ah, which is frequently used for the analysis and simulation of <100> 
silicon, are large. By combining published data for the low-field mobility with our data at higher 
fields, we determine parametrizations of the drift velocities of electrons and holes in <100> silicon 
for electric fields up to 50kV/cm and temperatures between 233 and 333 K. Their accuracy is 
about 5%, and they are available for simulating sensors built on high-ohmic <100> silicon. 
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6. Appendices 

6 .1. Appendix 1: Comparison of mobility parametrization for <111> silicon 

In this appendix we compare fits by the new KS parametrization (Eq. and the frequently used 
CT parametrization (Eq. to the tof measurements of Canali [15] for <111> silicon. Selected 
results are shown in Fig.j^ where the inverse mobility, l/fj.{E), is plotted versus electric field. 
The constant low field mobility is clearly visible, with a sudden transition to a linear increase of 
l/p. with electric field, which occurs around 1.8kV/cm for electrons and 2.5kV/cm for holes at 
300 K. The transition point decreases with decreasing temperature. From the lower part of Fig. [ 7 ] 
which shows the ratio of the fit results to the data, we observe systematic deviations between the 
data and both fits. However, the new KS parametrization shows somewhat smaller deviations, in 
particular at the transition between low- and high-held mobilities. 



Figure 7: Inverse mobility, 1/p,, of electrons (black) and holes (red) for <111> silicon at (a) 
245 K, and (b) 300 K. The data of Canali [H] (crosses) were Htted by the KS parametrization 
(Eq.|- solid lines) and the CT parametrization (Eq. — dashed lines). At the bottom the 
ratio of the Canali data and the values obtained from the ht of the KS (crosses) and the CT 
parametrization (diamonds) are shown. 

To quantify the difference between the hts and the data, we use the rdev (Eq.|^. The values of 
rdev are 0.22 (0.17) % for electrons, and 0.38 (0.28) % for holes at T = 245 (300) K for the KS ht. 
For the CT ht they are 0.40 (0.41) % for electrons and 0.47 (0.42) % for holes. We conclude that 
the KS parametrization provides a somewhat better description of Canali data for <111> silicon 
than the standard CT parametrization. 


6.2. Appendix 2: Current transients for the p-type sensor 

When performing the hts we were unable to achieve a good description of the current transients 
of the HPK p-type diode when we assumed an electric held, corresponding to a uniform p doping 
in the entire volume of the diode. However, adding a held-free region of 1.6 pm at the n+p side 
and one of 1 pm at the p’^'p side, resulted in a good description. To check the results we performed 
also measurements with the TCT setup at CERN-SSD Lab. It uses a wide-band ampliher and an 
Agilent scope with 2.5 GHz bandwidth and a sampling rate of 20GS/s. Using the Fast-Fourier- 
Transform method we successfully determined the transfer function, and performed hts to the 
measured current transients. Fig.shows the results for the HPK diode at 150 and 1000 V 
measured at 293 K. The data conhrm the conclusions from the Hamburg measurements that held- 
free regions close to the diode surfaces have to be introduced in the simulation of the n+p diode 
to obtain a satisfactory description. This is particularly the case for the front-side illumination 
shown in Fig. 8b where the signal is dominated by the drift of holes. Values in the range of 1 to 
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(a) 


(b) 


Figure 8: Comparison of the current transients of the p-type HPK diode measured with the CERN- 
SSD TCT setup using 675nm laser light at 293K (crosses) with the two simulations (lines). Solid 
lines: field-free regions of 1.6 pm at the n+p side and of 1pm at the p+p side. Dashed lines: 
without field-free regions. The results at 150 V are shown in red, and the ones at 1000 V in black, 
(a) Rear-side illumination with the laser light, and (b) front-side illumination. 


2 pm, which result in a good description by the simulation, appear reasonable. However, we do 
not understand why a field-free region is required for the n+p and not for the p^n diodes. 
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